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This work focuses on determining the latitudinal structure of ammonia vapor in Saturn's cloud layer near 
1.5 bars using the brightness temperature maps derived from the Cassini RADAR (Elachi et al. [2004], 
Space Sci. Rev. 115, 71-110) instrument, which works in a passive mode to measure thermal emission 
from Saturn at 2.2-cm wavelength. We perform an analysis of five brightness temperature maps that 
span epochs from 2005 to 2011, which are presented in a companion paper by Janssen et al. (Janssen, 
M.A., Ingersoll, A.P., Allison, M.D., Guilds, S., Laraia, A.L., Baines, K„ Edgington, S., Anderson, Y„ Kelleher, 
I<„ Oyafuso, F. [2013], Icarus, this issue). The brightness temperature maps are representative of the spa- 
tial distribution of ammonia vapor, since ammonia gas is the only effective opacity source in Saturn's 
atmosphere at 2.2-cm wavelength. Relatively high brightness temperatures indicate relatively low 
ammonia relative humidity (RH), and vice versa. We compare the observed brightness temperatures to 
brightness temperatures computed using the Juno atmospheric microwave radiative transfer (JAMRT) 
program which includes both the means to calculate a tropospheric atmosphere model for Saturn and 
the means to carry out radiative transfer calculations at microwave frequencies. The reference atmo- 
sphere to which we compare has a 3x solar deep mixing ratio of ammonia (we use 1.352 x 10 4 for 
the solar mixing ratio of ammonia vapor relative to H 2 ; see Atreya [2010], In: Galileo's Medicean Moons 
- Their Impact on 400 years of Discovery. Cambridge University Press, pp. 130-140 (Chapter 16)) and is 
fully saturated above its cloud base. The maps are comprised of residual brightness temperatures— 
observed brightness temperature minus the model brightness temperature of the saturated atmosphere. 

The most prominent feature throughout all five maps is the high brightness temperature of Saturn’s 
subtropical latitudes near ±9° (planetographic). These latitudes bracket the equator, which has some of 
the lowest brightness temperatures observed on the planet. The observed high brightness temperatures 
indicate that the atmosphere is sub-saturated, locally, with respect to fully saturated ammonia in the 
cloud region. Saturn’s northern hemisphere storm was also captured in the March 20, 201 1 map, and 
is very bright, reaching brightness temperatures of 166 K compared to 148 I< for the saturated atmo- 
sphere model. We find that both the subtropical bands and the 2010-2011 northern storm require very 
low ammonia RH below the ammonia cloud layer, which is located near 1.5 bars in the reference atmo- 
sphere, in order to achieve the high brightness temperatures observed. The disturbances in the southern 
hemisphere between -42° and -47° also require very low ammonia RH at levels below the ammonia 
cloud base. Aside from these local and regional anomalies, we find that Saturn's atmosphere has on aver- 
age 70 ± 1 5% ammonia relative humidity in the cloud region. We present three options to explain the high 
2.2-cm brightness temperatures. One is that the dryness, i.e., the low RH, is due to higher than average 
atmospheric temperatures with constant ammonia mixing ratios. The second is that the bright subtrop- 
ical bands represent dry zones created by a meridionally overturning circulation, much like the Hadley 
circulation on Earth. The last is that the drying in both the southern hemisphere storms and 2010- 
201 1 northern storm is an intrinsic property of convection in giant planet atmospheres. Some combina- 
tion of the latter two options is argued as the likely explanation. 

© 2013 Elsevier Inc. All rights reserved. 


1. Introduction 
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atmosphere for the past decade. Ideally, we would like to get a 
comprehensive picture of Saturn’s atmosphere that reconciles the 
general circulation, the cloud and haze distributions and composi- 
tions, the zonal wind profile, and the storm locations and dynam- 
ics. One major observational roadblock is that the stratospheric 
and upper tropospheric clouds and hazes on Saturn block our view 
of the atmosphere beneath them. 

The location and magnitude of the zonal jets at the cloud tops 
are well known from Voyager measurements (Sanchez-Lavega 
et al., 2000). The broad, strongly superrotating jet centered on 
the equator is a distinctive feature, with alternating eastward 
and westward jets to either side of the equator. Unlike Jupiter, con- 
vection on Saturn appears in both westward and eastward jets (Del 
Genio et al., 2009). Convective events on Saturn are intermittent, 
and the cause of the intermittency is uncertain. Saturn electrostatic 
discharges, or SEDs (Kaiser et al., 1983; Porco et al., 2005; Fischer 
et al., 2006, 2007), have been observed in convective storms and 
are indicative of lightning at depth. What causes these convective 
outbursts on Saturn, and how do they contribute to or maintain the 
general circulation? How does deep convection work on Saturn, 
and how does it fit together with the latitudinal belt-zone struc- 
ture of the giant planets? Answers to these questions have been 
difficult to obtain. The 2.2-cm observations analyzed in this work 
provide new data on the distribution of ammonia vapor in and be- 
neath the ammonia clouds, and will help diagnose the atmospheric 
dynamics at work inside the convective storms. 

The structure of Saturn’s clouds and hazes is still being studied, 
although the general features are understood. The equatorial zone 
is a region of constant high clouds and thick haze, whereas the 
midlatitudes (generally between ±20° and ±60°) are regions of 
smaller, more variable clouds (West et al., 2009). The vertical 
structure and composition of these clouds and hazes is not well 
known, but Cassini observations made by the ISS (imaging science 
subsystem), VIMS (Visual and Infrared Mapping Spectrometer) and 
CIRS (Composite Infrared Spectrometer) instruments are closing 
our knowledge gaps in these areas. Tied to the distribution of 
clouds and hazes is the distribution of tropospheric gases, for 
example ammonia and phosphine. How does the latitudinal distri- 
bution of clouds, hazes, and tropospheric gases coincide with Sat- 
urn’s belt-zone structure? Knowing the spatial distribution of 
these gases can help us determine the dynamical mechanisms that 
produce the spatial patterns themselves. For example, vertical mo- 
tion, caused by either convection or large-scale meridional over- 
turning, plays a key role in determining where clouds and hazes 
will or will not form. 

This work focuses on determining the latitudinal structure of 
ammonia vapor in Saturn’s ammonia cloud layer using the bright- 
ness temperature maps derived from the Cassini RADAR (Elachi 
et al., 2004) instrument, which works in a passive mode to measure 
thermal emission from Saturn at 2.2-cm wavelength. These maps 
are presented in a companion paper by Janssen et al. (2013, this is- 
sue), hereafter referred to as J13. The maps provide data on the 
spatial distribution of ammonia vapor in the pressure range 1-2 
bars, in the vicinity of the ammonia ice cloud. We believe these 
maps provide information about Saturn’s meridional circulation. 
The 2.2-cm data have better spatial resolution and sensitivity than 
any other microwave data on Saturn. The calibration of Cassini’s 
RADAR instrument, described in detail in Janssen et al. (2009) 
and J13, is accurate and was validated using both Saturn and more 
recent Titan observations as described inJ13. 

Section 2 describes the 2.2-cm observations and the radiative 
transfer model used in our analysis. The brightness temperature 
maps are described in Section 3. Section 4 compares the observa- 
tions to the output from the radiative transfer model. Discussion 
and implications for Saturn’s atmospheric dynamics are given in 
Section 5, and conclusions are given in Section 6. 


2. Observations and radiative transfer model 

Cassini’s RADAR radiometer was used to map Saturn during five 
equatorial periapsis passes occurring between 2005 and 2011. The 
maps were formed from continuous pole-to-pole scans taken 
through Saturn nadir during the periapsis passes, allowing the 
rotation of Saturn to sweep the scan westward in longitude. The 
observations and mapping are described in detail in J13 along with 
the calibration and error analysis. We refer the reader to Section 2 
of J13 for a description of the observations and observational ap- 
proach, and to Section 3.2 of J13 for a description of the map-gen- 
erating process. 

The reference model used to calculate the residual brightness 
temperature maps is also described in detail in Section 3.1 of J13. 
The model and radiative transfer calculations were made using 
the Juno atmospheric microwave radiative transfer (JAMRT, Jans- 
sen et al., 2005, in preparation) program, which is in development 
for the Juno Microwave Radiometer (MWR) experiment on Jupiter. 
To match the RADAR observations, radiative transfer calculations 
are carried out at 2.2-cm wavelength (13.78 GHz), and brightness 
temperatures are output for each observation. This model builds 
an atmosphere with user-prescribed physical parameters, such as 
the vertical mixing ratio profiles of ammonia, phosphine and 
water. Temperature and pressure profiles are calculated assuming 
hydrostatic equilibrium using both wet and dry adiabats. The ref- 
erence model assumes a moist adiabatic temperature profile with 
100% relative humidity (RH), with a dry adiabatic profile below 
cloud base, such that the temperature is monotonically decreasing 
from the bottom to the top layer of the model atmosphere. The 
adiabats include the contributions from the NH 4 SH and H 2 0 
clouds, although the weighting function drops to essentially zero 
before we reach the water cloud at great depth. A temperature of 
134.8 K (Lindal et al., 1985) is specified at a pressure of 1 bar, 
and the model temperature profile is slaved to this reference value. 
We varied this value in order to test the sensitivity of the 2.2-cm 
brightness temperature to variations in the 1-bar temperature, 
and found the brightness temperature to be only minimally sensi- 
tive to this reference value (see Section 5.1). The topmost level of 
the model is the level at which the temperature reaches 1 10 K, 
which is 560 mb for the (134.8 K, 1 bar) reference point. The model 
assumes a completely transparent atmosphere above 110 K and 
therefore ignores this region of the atmosphere. The deepest level 
of the model atmosphere is 1000 bars, which is well below the 
pressure level sensed by the 2.2-cm observations, and the vertical 
layers are 100 m thick. The model also includes the emission angle 
dependence (limb darkening) of the brightness temperature. 

Table 1 gives the atmospheric constituents and their respective 
abundances in the model atmosphere, including the values used 
for the solar abundances. H 2 0, NH 3 , PH 3 , and H 2 S are the condens- 
able gases (Atreya, 2010). H 2 S reacts with NH 3 to form an NH 4 SH 
cloud with a base around 5 bars. An ammonia ice cloud forms 
above this, with a base around 1.5 bars. The water cloud is deeper 
(base ~10 bars) and out of the sensitivity range of the 2.2-cm 

Table 1 

Abundances of atmospheric constituents in the JAMRT program. Solar and enrichment 
values are from Atreya (2010), who calculated solar abundances from the photo- 
spheric values of Grevesse et al. (2005). 


Constituent 

Solar abundance (relative to H 2 ) 

Enrichment relative to solar 

He 

0.195 

0.6955 

ch 4 

5.50 x 10~ 4 

9.4 

h 2 o 

1.026 x 10- 3 

3.0 

nh 3 

1.352 x 10~ 4 

3.0 

h 2 s 

3.10 x 10~ 5 

5.0 

Ar 

7.24 x 10~ 6 

1.0 

ph 3 

5.14 x 10~ 7 

7.5 
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weighting function. In the model, the presence of the ammonia ice 
cloud particles does not affect the 2.2-cm brightness temperature 
significantly, although the depletion of ammonia vapor by the for- 
mation of the clouds does. Using Cassini Visible and Infrared Imag- 
ing Spectrometer (VIMS) data, Fletcher et al. (2011a) find evidence 
of a compact cloud deck in the 2. 5-2. 8 bar region. They point out 
that this is “deeper than the predicted condensation altitudes for 
pure NH 3 clouds (1.47-1.81 bars) and higher than the predicted 
condensation altitudes for the NH 4 SH cloud (4.56-5.72 bars).” 
They conclude, “The VIMS 2. 5-2. 8 cloud cannot be identified 
unambiguously using the present data set.” Resolving this discrep- 
ancy is beyond the scope of this paper, so we use an equilibrium 
condensation model similar to that referred by Fletcher et al. 
(2011a), even though it does not entirely fit their interpretation 
of the VIMS data. For a more in-depth description of the JAMRT 
program, see Section 3.1 of J13. 

Because we only have data at one wavelength, there arises the 
unavoidable ambiguity of whether temperature or ammonia is the 
cause of the variations in brightness temperature. Our analysis as- 
sumes that atmospheric temperature is constant with latitude in 
the sensitivity range of 1-2 bars, and that variations in the ammo- 
nia mixing ratio cause brightness temperature variations. From a 
data-fitting point of view, one could also perform this analysis 
assuming that the mixing ratio of ammonia is constant with lati- 
tude and that atmospheric temperature fluctuations cause bright- 
ness temperature variations. In general, brightness variations can 
be due to both fluctuations in ammonia concentration and atmo- 
spheric temperature from latitude to latitude, and these are not 
easily separated. We think that ammonia dominating the bright- 
ness temperature variations at 2-cm wavelength is the right 
choice, based on both the fact that emission from ammonia within 
the cloud region is strongly buffered against temperature varia- 
tions, and also the extreme sensitivity of condensation/evaporation 
to vertical motions, as exhibited in the Earth’s tropics. Large-scale 
subsidence, for example, will cause a “drying” of the atmosphere at 
certain latitudes, and could potentially produce a thermal emission 
pattern like the one we observe on Saturn. Our interpretation of 
the brightness temperature variations as variations in ammonia 
abundance is consistent with Grossman et al. (1989), who analyzed 
thermal emission from Saturn at 2- and 6-cm wavelength. They ar- 
gue that large temperature deviations (on the order of 8 K in their 
paper) would be difficult to sustain in the presence of convection. 

There are two basic ways to increase the brightness tempera- 
ture Tb. Either increase the atmospheric temperature T with con- 
stant ammonia mixing ratio, or else hold T constant and lower 
the ammonia mixing ratio. In both cases the relative humidity 
(RH) goes down. If RH stays constant, an increase (decrease) in 
temperature is offset by an increase (decrease) in ammonia abun- 
dance due to the Clausius-CIapeyron relation, and the brightness 
temperature stays the same. Thus Tb is measuring RH in the sense 
that high RH gives low Tb and vice versa. When we refer to “drying 
out” of the atmosphere, we are referring to low RH acting to pro- 
duce high Tb. 

We varied two parameters of the model to produce departures 
from the reference model: the enrichment factor of ammonia rela- 
tive to the solar abundance (EF), and the ammonia depletion factor 
(DF). The enrichment factor EF is defined as the deep mixing ratio 
of ammonia vapor, expressed in terms of the solar abundance of 
ammonia (Table 1). In the model ammonia is uniformly mixed be- 
low the level where it reacts with H 2 S to form an ammonium 
hydrosulfide cloud. It is partially depleted from that level up to 
the ammonia condensation level, where the mixing ratio of ammo- 
nia falls off according to the saturation vapor pressure dependence 
on temperature (the Clausius-CIapeyron relation). The deep abun- 
dance of ammonia is not precisely known on Saturn, but is thought 
to be in the range 2-4x solar (Atreya, 2010). We vary EF from 2 to 


8 in this work (corresponding to volume mixing ratios of 2.3- 
9.4 x 1CT 4 ). This range encompasses previous estimates for the 
deep abundance of gaseous ammonia, for example 4-6 x 10~ 4 re- 
ported by Briggs and Sackett (1989) and 5 x 1CT 4 reported by de 
Pater and Massie (1985) from VLA measurements. 

The depletion factor DF X is chosen to allow an additional deple- 
tion of ammonia above some level X. Beginning with an ammonia 
mixing ratio distribution determined for an atmosphere with a gi- 
ven EF, DF X is simply a scale factor between 0 and 1 that multiplies 
the vertical distribution of ammonia above level X. It is intended 
purely as a first-order parameter to investigate ammonia deple- 
tions likely to occur in more realistic dynamical atmospheres. In 
applying DF X we ignore any perturbations implied for other com- 
ponents of the model atmosphere such as the cloud base height 
or the temperature profile. For example, in the cloud layer, DF cb 
would correspond to the relative humidity (RH) of ammonia, 
where we choose the cloud base as our level X. In the following 
we choose DF cb and a second choice DF 5bar , with level X as the 5- 
bar pressure level, to effectively bracket the cases we will study 
for ammonia depletion. We will show that both parameters are 
needed to explain the observed 2-cm brightness temperatures. 

Figs. 1-3 demonstrate the effects of varying EF and DF in the 
model atmosphere. Fig. 1 shows two ammonia profiles, given by 
the heavy solid (EF = 8) and dashed (EF = 3) lines, both with 
DF = 1 (at all altitudes). The ammonia mixing ratios are less than 
3 and 8x solar at the 2.5-bar limit of Fig. 1 because the NH 4 SH 
cloud has a base around 5 bars and depletes the ammonia above. 
The light solid and dashed lines are the weighting functions that 
correspond to the 8x solar and 3x solar ammonia profiles, respec- 
tively. The weighting function for 3x solar extends deeper because 
there is less ammonia at the 1.5-1. 8 bar level to block the radiation 
from below. The dotted curves are the temperature profiles for the 
two atmospheres, which differ only very slightly at the 2.5 bar le- 
vel (with the 3x solar case being slightly warmer than the 8x solar 
case). The calculated brightness temperatures for these two 
models are 148.0 K for 3x solar and 147.9 K for 8x solar ammonia. 
Increasing EF further has very little effect on the brightness 
temperature. 


Temperature (K) 

120 140 160 180 



Fig. 1 . Two vertical profiles of ammonia vapor with varying EF (bold lines). The 
solid and dashed bold lines are profiles with EF = 8 and 3, respectively, and the solid 
and dashed lines are their respective weighting functions. Their temperature 
profiles (top x-axis) are given by the dotted lines, which are almost identical in this 
pressure range, except that the 3 x solar case is slightly warmer than the 8 x solar 
case at 2.5 bars. The model brightness temperatures for the 3 and 8x solar cases are 
148.0 1< and 147.9 K, respectively. 
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Temperature (K) 



EF (NH 3 /solar) 


Fig. 2. Same as Fig. 1, but for constant EF = 3 and varying DF cb (or ammonia RH in 
the cloud layer). In this case the two model atmospheres have the same 
temperature profile. The case with DF cb = 1 is identical to the EF = 3 case in Fig. 1 
(bold dashed line). There is a discontinuity in the ammonia mixing ratio at the cloud 
base for DF cb = 0.2 (bold solid line) due to the way the DF parameter functions in the 
model. The model brightness temperatures are 148.0 K and 154.0 1< for the DF cb = 1 
and 0.2, respectively. 


Temperature (K) 

120 140 160 180 200 220 240 



EF (NH 3 /solar) 


Fig. 3. Same as Figs. 1 and 2, but for constant EF = 3 and varying DF 5bar . Again, the 
DF 5bar = 1 case is the same as EF = 3 in Fig. 1 and DF cb = 1 in Fig. 2. Note the change 
in scale of the y-axis. Like Fig. 2, there is a discontinuity in the ammonia mixing 
ratio at the level at which DF is applied (bold solid line), which is 5 bars for this case. 
The ammonia continues to be depleted between 4 and 5 bars because the formation 
of the NH 4 SH cloud at 5 bars causes ammonia depletion there. The model 
brightness temperatures are 148.0 K and 161.2 K for DF 5bar = 1 and 0.2, respectively. 


Fig. 2 is like Fig. 1 but with EF = 3 held constant and DF cb varied. 
The heavy solid and dashed lines are for atmospheric profiles of 
DF cb = 0.2 (RH = 20%) and 1 (RH = 100%), respectively. The dotted 
line is the temperature profile of both atmospheres. The calculated 
brightness temperatures for DF cb = 1 and 0.2 are 148.0 K and 
154.0 K respectively. 

Fig. 3 is like Fig. 2, except DF 5bar is varied. Note the difference in 
scale of the y-axis between Fig. 3 and the previous two figures. In 
addition to the imposed depletion above the 5 bar level, the NH 4 SH 
cloud also depletes ammonia above the 5 bar level, which is why 


there is still some ammonia depletion above 5 bars for DF 5bar = 1. 
The calculated brightness temperatures for DF 5bar = 1 and 0.2 are 
148.0 K and 161.2 K, respectively. 

The explanation for the differing brightness temperatures in 
each model atmosphere is as follows: If there is relatively less 
ammonia in a given atmospheric column (from the top down), 
then the 2.2-cm weighting function will have contributions from 
higher pressures. Because the model temperature falls off adiabat- 
ically with height at all levels, the 2.2-cm brightness temperature 
will be higher when there is less ammonia in the column, since 
we are probing a lower altitude in the atmosphere where the tem- 
perature is warmer. The opposite is true if there is relatively more 
ammonia in the column; namely, the resulting brightness temper- 
ature will be lower. 

We chose to vary DF down to a level of 5 bars, which sparks the 
question: how deep is it necessary to deplete ammonia in order to 
achieve the observed brightness temperatures? Fig. 4 helps answer 
this question by displaying DF X as a function of X, where X is the 
pressure level to which we deplete ammonia (EF = 3 for all of these 
calculations). Except for the smallest values of DF X , the curves are 
flat forX > 2 bars. Depleting down to 5 bars is equivalent to deplet- 
ing down to 2 bars, which means that the 2.2-cm weighting func- 
tion is very small below 2 bars in this parameter regime (EF = 3, 
0.1 c DF X Cl), and the brightness temperature is not very sensi- 
tive to depletion below the 2 bar level. For values of DF X less than 
0.1, the brightness temperature is very sensitive to the depth of 
depletion because the weighting function peaks at X bars for small 
DF X . A similar plot was made for EF = 6, and the difference was that 
the brightness temperature of the 6x solar case was less than the 
3x solar case by 1-7 1< depending on the value of DF X . The lower 
brightness temperatures for EF = 6 are to be expected because 
when EF is larger, there is less emission from the deeper levels 
and the upper levels must be depleted a little more to get the same 
brightness temperature. 

3. Maps 

After using the radiative transfer model to calculate the 
brightness temperature as a function of emission angle, the 
residual brightness temperatures were calculated relative to a 
saturated atmosphere (RH = 100% above the condensation level) 



Fig. 4. Contours of model brightness temperature (in Kelvin) as a function of the 
depletion factor DF X and X (the depth of depletion) for constant enrichment factor 
of ammonia relative to the solar abundance EF = 3. DF X is a scalar parameter 
between 0 and 1 that multiplies the vertical distribution of ammonia above some 
pressure level X. Thus, DF X = 1 is no depletion and DF x = 0 is 100% depletion of 
ammonia vapor above level X. 
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model with the constituent enrichments given in Table 1. All future 
mention of residual brightness temperatures refers to the residuals 
from this 3x solar ammonia reference model, which was chosen 
based on the ability of the EF = 3 reference model to span the ob- 
served brightness temperatures better than the models with high- 
er or lower values of EF (Section 4, Figs. 1 1 and 12). Fig. 5 is a map 
from December 9, 2009 of the 2.2-cm residual brightness temper- 
ature, created by the procedure explained in Section 3.2 of J13. 
There are four other maps analyzed in this work but not shown 
here: September 23, 2005, October 13, 2009, July 25, 2010 and 
March 20, 2011. They are presented in Fig. 9 of J13. 

3.3. Obstruction by the rings 

An issue persistent throughout all the maps is that the equator 
is obstructed by the rings (black region along equator in Fig. 5). The 
rings are optically thick scatterers with very little intrinsic thermal 
emission, and hence lower the measured brightness temperature 
in the northern (southern) hemisphere when the spacecraft is be- 
low (above) the ring plane. When the spacecraft is in the ring 
plane, the ring inclination angle is exactly 0° and the rings disap- 
pear from view, allowing full view of the equator. Fig. 6 demon- 
strates the effect of the ring obstruction on brightness 
temperature for the July 2010 map, the map for which the effect 
of the rings is most prominent (see Fig. 9 of J13). Each line is a dif- 
ferent meridional (pole-to-pole) scan near the spacecraft’s ring 
plane crossing (RPC), labeled with the ring inclination angle. For 
an example of a RPC scan see Fig. 5 near longitude 18°W, where 
the ring blockage at the equator goes to zero. During the July 
2010 observation period, Cassini was moving rapidly across the 
ring plane. Thus the effect of the rings on the brightness tempera- 
ture is large, 50% between two scans taken only minutes apart (e.g. 
from scan -0.019° to -0.377°). Fig. 6 demonstrates that ring incli- 
nation angles as small as 0.1° have large effects on the observed 
equatorial brightness temperature. 

Cassini crossed the ring plane in four of the five maps. We test 
our ability to remove the ring effect with a model that assumes an 
isothermal (150 K) brightness temperature for the atmosphere and 
a ring brightness temperature of 25 K. With this model we 
synthesize the individual brightness scans taking into account 
the actual geometry of the spacecraft and the instrument. The right 
panel of Fig. 7 shows the results from this model for the four scans 



Planetographic Latitude (degrees) 

Fig. 6. Five individual scans around ring plane crossing (RPC) for the July 2010 map. 
Scans are labeled with their ring inclination angles as viewed by Cassini. The first 
scan, labeled 0.318°, was made while the spacecraft was in the southern 
hemisphere, therefore the ring blockage occurred in the northern hemisphere. As 
Cassini approached the ring plane, the effect of the rings became quite small. 
Because the spacecraft was moving fairly quickly across the ring plane in this map, 
the rings had a very large effect from one scan to the next. 


closest to RPC for the December 2009 map, which has the smallest 
ring inclination angle while observing the equator. The labels are 
the same as in Fig. 6, but the scales along the y-axis are different. 
The left panel shows the observed brightness temperatures for 
the same scans. The model does not predict any dip in brightness 
temperature at the equator for the RPC scan (0.001° ring 
inclination angle), demonstrating that we are observing the true 
equatorial brightness temperature in this scan. In the left panel, 
the 0.001° scan is continuously flat through the equator (within 
±2°), which does not occur for any other RPC scan except for the 



360 330 300 270 240 210 180 150 120 90 60 30 0 

West Longitude 


Tr (K> 1 

0 8 16 

Fig. 5. 2.2 -cm residual brightness temperature (in Kelvin) map of Saturn from December 9, 2009. The residual temperature is calculated by subtracting the brightness 
temperature from a fully-saturated reference model with 3 x solar ammonia mixing ratio (Section 2) from the observed brightness temperature. The black band at the equator 
is due to the cold rings obstructing the atmosphere. Section 3.2 ofJ13 offers a detailed explanation of how the brightness temperature maps were generated. 
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Planetographic Latitude (degrees) Planetographic Latitude (degrees) 

Fig. 7. Left panel: Same as Fig. 6 for the December 2009 map. The RPC scan for this map has the smallest ring inclination angle of all the maps (while viewing the equator), 
and is flat across the equator. This is the best view of the equatorial brightness temperature that we have of all five maps. Right panel: Same labeling, but for a simple beam 
convolution model that takes into account the Cassini-Saturn geometry and includes the A and B rings only. It assumes an isothermal atmosphere of 1 50 K, and takes the 
microwave brightness of both the A and B rings to be 25 K. According to this model, we actually see the equator with no ring blockage for the 0.001° scan. 



Fig. 8. Residual brightness temperatures for RPC scans of the four maps and their 
average. The ring inclination angle for each scan is in parentheses next to the date. 
Because of the low resolution of the September 2005 scan, the bright bands off the 
equator affect the equatorial brightness temperature, causing it to be 0.7 K too high. 
With this in mind, the four scans are within ±1.5 K. The December 2009 RPC scan is 
flat across the equator with a very small ring inclination angle. Thus it provides the 
true equatorial brightness temperature. 
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Fig. 9. Mean residual brightness temperature (open circle and + signs) and its mean 
standard deviation as a function of planetographic latitude (solid line) for all five 
maps, excluding the northern storm and the latitudes near the equator where the 
rings block the view of the atmosphere. Observations were sorted into latitude bins 
0.4° wide since the latitudes are unevenly spaced. From -4° to +4° the single 
December 2009 RPC scan is used (+ signs), because it is the best view of the 
equatorial brightness temperature that we have of all the maps (Section 3.1). 
Standard deviations were calculated at each latitude for each map and then 
averaged. The average is weighted by the number of observations at a given latitude 
and date and the sum is over all five observation dates. 


September 2005 scan, which is at very low resolution compared 
with the other three maps (Fig. 8). Thus we take the December 
RPC scan as the RPC scan that provides the true equatorial bright- 
ness temperature. The increase in brightness temperature at lati- 
tudes greater than 3° in both hemispheres is not in the model 
but is a real property of Saturn’s atmosphere. 

Fig. 8 displays residual brightness temperatures versus latitude 
for all four RPC scans. The scale along the y-axis is expanded 
relative to that in Figs. 6 and 7. The ring inclination angles for each 
scan are displayed in parentheses in the legend. The segment of the 


September 2005 map where the RPC occurred is at very low reso- 
lution (Fig. 9,J13), and is therefore affected by the two bright bands 
at ±9° with a contribution of about +0.7 K. Applying this correction 
to the September 2005 scan brings the residual brightness temper- 
ature down to 0.5 K at the equator. The observed residuals are con- 
sistent with ~1 K variability in the equatorial region. We chose the 
December 2009 RPC scan to represent the equatorial brightness 
temperature for all five maps in Fig. 9 because it has the best 
geometry (lowest inclination at RPC) and produces the best picture 
of the equatorial brightness temperature that we have. 
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3.2. General features 

With the exception of the March 2011 map, which has the 
2010-2011 northern storm (Figs. 9 and 12, J1 3) or great white spot 
(Fischer et al„ 2011; Sanchez-Lavega et al., 2011), all of the maps 
share the same general characteristics. In what follows, all lati- 
tudes are planetographic unless explicitly stated. In Fig. 5, the 
equatorial region, within 10° of the equator, is texturally anoma- 
lous compared with the rest of the map, even when excluding 
the effect of the rings. There is non-uniform high brightness near 
-9° and 9° that generally decreases towards the equator, which 
is obstructed by the rings for all scans except for the one during 
RPC. The brightness temperature variations on the maps are quite 
large, with variations of more than 10K, and we investigate the 
causes of these variations. There are some structures in the south- 
ern hemisphere band between -42° and -47° (e.g. at 340°W and 
75°W in Fig. 5), which is just south of the westward jet at -42° 
(-35° planetocentric) and has been the site of many lightning 
observations (Dyudina et al., 2007; Fischer et al., 2011). There are 
also two narrow bright bands at -33° and -37°, and a broad dark 
band from -15° to -30°. 

The highest brightness temperature in the northern storm is 
165.7 K (18.9 K residual brightness temperature). The highest 
brightness temperature in all five maps is 167.1 K (19.3 K residual 
brightness temperature), which occurs in the subtropical latitudes 
in the October 2009 map. 

Fig. 9 shows the zonally averaged residual brightness tempera- 
ture (open circles and + signs) and its standard deviation (solid 
line) as a function of latitude. Outside of ±4° latitude, the bright- 
ness temperatures are averaged over longitude and over all five 
maps (open circles), excluding the 2010-2011 northern storm 
and the latitudes near the equator where the rings block the view 
of the atmosphere. Within ±4° of the equator, a single RPC scan 
from the December 2009 map is used (+ symbols), as explained 
in Section 3.1. The standard deviation was computed with respect 
to longitude for each latitude outside ±4°, and then the weighted 
standard deviation for all five maps was calculated. The weights 
used in this calculation are the number of observations at a given 
latitude and date and the sum is over all five observation dates. 
The globally averaged residual brightness temperature from 
Fig. 9 is 1.7 ± 1.1 K, where the 1.1 K is real variability of the longi- 
tudinally-averaged brightnesses in Saturn’s atmosphere. 

The residual brightness temperatures shown in Fig. 9 are posi- 
tive at every latitude. One important implication of this observa- 
tion is that the atmosphere is always ammonia-depleted with 
respect to a fully-saturated 3x solar ammonia model (Table 1). A 
striking feature of Fig. 9 is the two relatively bright bands near 
±9°, the subtropics of Saturn, with residuals of 3.8 l< and 5.8 K in 
the southern and northern hemispheres, respectively (correspond- 
ing to brightness temperatures of 151.7 K and 153.7 K). These two 
bright latitudes surround a relatively low residual brightness tem- 
perature of 0.1 K (corresponding to a brightness temperature of 
148.1 K) at the equator. At the equator the atmosphere is close to 
being saturated with ammonia. The subtropical bands are accom- 
panied by elevated standard deviation, indicating that there is 
some structure in these regions. Another feature is the pair of 
bright bands around -36° and -34° with averaged residuals of 
2.5 and 3.6 K respectively (corresponding to brightness tempera- 
tures of 149.2 and 150.5 l<), which correspond to the two narrow 
bands in the southern hemisphere in Fig. 5. There are two regions 
with high standard deviation in the southern hemisphere, one near 
-33° (-28° planetocentric) and the other near -43° (-37° plane- 
tocentric). The latter is the latitude of storm alley and the site of 
the southern hemisphere lightning (Dyudina et al., 2007). The rel- 
atively high standard deviation in storm alley corresponds to the 
bright dots seen there (Fig. 5, longitude 345°W), which are likely 


to be holes in the ammonia layer associated with the holes in the 
clouds described by Dyudina et al. (2007). Dyudina et al. (in prep- 
aration) investigates the structure of these southern hemisphere 
storms and presents the lightning observations from both the 
southern storms and the 2010-2011 northern storm. The northern 
latitudes have relatively constant brightness temperatures, with 
fluctuations on the order of 1 K from latitude to latitude. The 
southern hemisphere has larger brightness temperature gradients 
than the northern hemisphere, for example almost a 4 K increase 
from -25° to -35°. For the northern storm (not included in 
Fig. 9), the standard deviation in the latitude band between 20° 
and 50° is much larger, reaching a peak value of 6 l< at 40°. 

When looking at the relatively bright spots within the bright 
subtropical bands in all five maps, it is natural to ask if there is 
any periodic structure in these regions. Jupiter, for example, had 
the equatorial plumes in the northern hemisphere with longitudinal 
wavenumber between 1 1 and 1 3 at the time of the Voyager encoun- 
ters (Allison, 1990), and between 8 and 12 determined more re- 
cently (Arregi et al., 2006). To determine whether there is periodic 
structure in the subtropical bands of Saturn, we calculated the auto- 
correlation of the brightness temperature with respect to longitude 
for latitude bands between 6° and 1 0° in both hemispheres. We used 
binning to remedy the problems that arise due to unevenly spaced 
observation points and large gaps in the maps. For a given latitude, 
every pair of points was placed into a bin 3° wide based on the lag 
between the points. Thus the 0° lag bin has pairs of points with lags 
from 0° to 3°. The next bin is 3-6°, and so on. We did this for four 
maps, excluding the July 2010 map because the rings obstruct the 
majority of the subtropical bands. The resulting averaged correla- 
tion coefficients are plotted versus longitudinal lag in Fig. 10. The 
correlations in the southern hemisphere tend to fall off more rapidly 
than in the northern hemisphere, with the exception of the Decem- 
ber 2009 map. This indicates that the bright spots in the northern 
hemisphere subtropical band have a greater longitudinal span than 
those in the southern hemisphere. Four of the eight panels in Fig. 10 
reveal wave-like features with longitudinal periods ranging from 
20° to 45° (zonal wavenumbers 18-8, respectively) for the 6-10° 
latitude range. The period varies from year to year, and there is no 
indication of a recurring dominant period. We did not find any 
wave-like features in other latitude ranges. 

4. Comparison with radiative transfer model 

Ammonia vapor is the only effective source of opacity in Sat- 
urn’s atmosphere at 2.2-cm wavelength, as shown by Fig. 3 of 
J13 that plots absorption coefficient versus height for the relevant 
constituents. Ammonia is by far the dominant absorber in our sen- 
sitivity range. Thus the 2.2-cm brightness temperature maps yield 
information about the ammonia vapor distribution. To interpret 
the maps we used the JAMRT program described in Section 2 to 
calculate 2.2-cm brightness temperatures based on different verti- 
cal profiles of ammonia. We varied two parameters, EF and DF, also 
described in Section 2. 

Figs. 1 and 2 suggest that we cannot get brightness temperature 
variations larger than 11 K by varying only EF (ammonia deep 
abundance) or DF cb (ammonia RH in the cloud layer), since these 
two parameters only account for brightness temperature varia- 
tions of up to 11 K. Fig. 11 demonstrates this by displaying the 
model brightness temperature (in Kelvin) as a function of DF cb 
(or RH) and EF. All the calculations done in Figs. 11 and 12 were 
done at an emission angle of 0°. This is important to note because 
the model brightness temperatures decrease with increasing emis- 
sion angle due to limb darkening. For example, a brightness tem- 
perature of 160 K at the equator would have a residual of 12 K, 
while a brightness temperature of 160 K at 40° (e.g., the 2010- 
2011 northern storm) would have a residual of 13.5 K. 
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Fig. 10. Autocorrelations of brightness temperature with respect to longitude versus longitudinal lag. The panels display averages of the autocorrelations in the latitude 
bands from 6° to 10° in the northern hemisphere (left column) and the southern hemisphere (right column) for each map, excluding the July 2010 map due to the large 
amount of ring blockage in the subtropics. 


In Fig. 11, as expected, the brightness temperatures increase 
with decreasing DF cb (decreasing RH), since depletion of ammonia 
pushes the weighting function deeper into the atmosphere where 
it is warmer. However, as mentioned above, the largest tempera- 
ture contrast that can be obtained in this parameter regime is 
11 K (148-159 K, Fig. 11), which indicates that ammonia vapor 
must be depleted to levels deeper than the ammonia cloud in order 
to achieve the highest brightness temperatures seen in the maps 
(160+ I<), provided that DF cb > 0.1. If DF cb < 0.1, then the ammonia 
concentration in the cloud region is almost zero, which is unlikely 
given ammonia concentrations reported by Fletcher et al. (2011a) 
for the 1-4 bar region. The idea that the formation of the NH 4 SH 
cloud around 5 bars is the ammonia-depletion mechanism beneath 
the ammonia-ice cloud has been discussed by previous authors 
(e.g. Briggs and Sackett, 1989). Since our model includes the forma- 
tion of this cloud (with H 2 S enriched by 5x solar, Table 1), ammo- 
nia must be depleted even more than just by the NH 4 SH cloud 
formation in order to agree with observations, unless H 2 S is in real- 
ity enriched by more than 5x solar on Saturn. 


There is a trend of increasing brightness with increasing EF at 
low DF cb in Fig. 11. This is counter-intuitive, because increasing 
EF means more ammonia, but it is an artifact of the model and 
has an explanation. The DF cb parameter depletes ammonia down 
to the ammonia cloud base, and for very low DF cb that is where 
the weighting function peaks. As EF increases, the ammonia cloud 
base moves to deeper (and warmer) levels (Fig. 1), bringing the 
weighting function with it. The effect of lowering the altitude of 
the weighting function outweighs the effect of adding ammonia 
at deeper levels, and brightness temperature increases with 
increasing EF. 

Fig. 12 is the same as Fig. 11, but for DF 5bar . Note that the con- 
tour interval has been increased from 1 to 4 K. As expected, the 
brightness temperature increases with decreasing DF 5bar and 
decreasing EF. By depleting the ammonia down to this deeper level 
(5 bars), we are able to achieve the 166-167 K brightness temper- 
atures seen in the maps. This requires either DF cb <0.1 (Fig. 11, 
lower right portion) or a combination of DF 5bar <0.3 and EF<4 
(Fig. 12, lower left corner). Figs. 11 and 12 are for rays propagating 
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Fig. 11. Contours of model brightness temperature (in Kelvin) as a function of the 
ammonia depletion factor above the cloud base height, DF C b (or RH), and the 
enrichment factor of ammonia relative to the solar abundance, EF. The warming 
trend at high EF is related to the way the cloud base changes with this parameter 
(see text for a more detailed explanation). 



Fig. 12. Same as Fig. 11, but for the depletion factor above 5 bars, DF 5ba r. Note the 
contour interval has been increased to 4 K. 

vertically. Taking the limb darkening effect at 40° latitude into ac- 
count means that brightness temperatures in the storm correspond 
to DF 5bar = 0.1 for EF = 3, which is a very large depletion (90%). The 
subtropical bright bands can be explained by either very low DF cb 
or mid-range values of DF 5bar . 

From our results it seems that ammonia lies near the 3x solar 
range, because at high EF values (EF > 4) it becomes impossible 
to achieve the highest observed brightness temperature (167 K) 
unless DF 5bar < 0. 1 ( Fig. 12). Such low values of DF are incompatible 
with other observations of ammonia (i.e. Fletcher et al., 2011a). 
Also, at low EF values (EF<3) it becomes impossible to achieve 
the lowest observed brightness temperature (148 K at the equator) 
unless the atmosphere there is supersaturated. This is because 
saturation occurs at T>1481< for EF<3 (Figs. 11 and 12, top left 
corner). According to this analysis, Saturn’s atmosphere lies in 
the 3-4 x solar ammonia range, with fairly large depletion of 
ammonia extending below the cloud base in some regions. 

Since the highest brightness temperature seen in all five maps is 
near 167 I< (in the subtropical latitudes of the October 2009 map), 


depleting down to 2 bars could explain all the brightness temper- 
atures we see (Fig. 4). In this case DF 2bar would have to be close to 
0, which corresponds to no ammonia above 2 bars. This is not a 
likely scenario, so we focus on DF 5bar for the remainder of the pa- 
per, keeping in mind that DF 5bal - is the same as DF 3bar or DF 4bal - as 
long as DF ^ 0.2 (Fig. 4). 

5. Discussion 

The deep abundance of heavy elements such as nitrogen and 
carbon is not well known on Saturn. Estimates range from 2 to 
4x solar (Atreya, 2010) for ammonia to 9-10x solar (Fletcher 
et al., 2012) for carbon. The solar values used in these estimates 
are those given in Table 1, and they are given with respect to H 2 . 
Because we only have data at one wavelength, we cannot constrain 
EF and DF separately, we can only comment on possible combina- 
tions of the two parameters that give brightness temperatures con- 
sistent with the 2.2-cm data. However, as presented above, EF 
must lie in the 3-4 x solar range in order to achieve the highest 
and lowest brightness temperatures observed at 2.2-cm. For the 
following discussion, we assume that EF = 3 and then we comment 
on the values of DF that would yield the observed brightness tem- 
peratures. Note that EF = 3, which corresponds to a deep volume 
mixing ratio of 3.6 x 1CT 4 , is fairly consistent with values reported 
by de Pater and Massie (1985) and Briggs and Sackett (1989), 
whose estimates are for pressure levels greater than 3 bars and 
the 2-bar level, respectively, and Fletcher et al. (2011a) for the 1- 
4 bar range at the equator. Our estimate is also consistent with 
the estimate of 1.2 x 10~ 4 obtained from VLA measurements by 
Grossman et al. (1989) for the condensation altitude. 

5.3. All maps 

Using the quantitative analysis from the radiative transfer mod- 
el for 3x solar ammonia, we were able to convert the residual 
brightness temperatures to DF cb values, or values of ammonia RH 
in the ammonia cloud layer. Fig. 13a shows a map of the ammonia 
RH for the March 201 1 observation date, with the northern storm 
being the most prominent feature. Black regions would be regions 
where the atmosphere is supersaturated with respect to ammonia 
in the cloud layer (we do not see any). In this figure black corre- 
sponds to where the atmosphere looks very cold due to the ring 
obstruction. Blue regions are regions that require depletion of 
ammonia below the clouds in order to achieve those high bright- 
ness temperatures. 

Fig. 9 shows that the global average residual brightness temper- 
ature relative to the saturated model is 1.7 ± 1.1 K, where the 
±1.1 K is the real variability of longitudinally-averaged brightness- 
es in the atmosphere. Figs. 4, 11 and 12 give DF for various bright- 
ness temperature values, and the curves are almost flat for X > 2 
bars (Fig. 4). In Fig. 1 1, each DF value is the RH at the corresponding 
brightness temperature. Taking brightness temperature Tb = 148 K 
as the saturated case and Tb = 149.7 ± 1.1 I< as the global average 
(for EF = 3), we find that the average RH is 70 ± 15% in the cloud 
layer. Note that this estimate does not include the 2010-2011 
northern storm, and that the ±15% comes directly from the 1.1 K 
variability in Fig. 9. For comparison, de Pater et al. (2001 ) find that 
the disk-averaged relative humidity of ammonia in Jupiter’s atmo- 
sphere is of the order of 10% at pressures less than 0.55 bars. 

Saturn’s atmosphere also has local and regional features with 
RH < 0 in the cloud layer. These regions are shown in green-blue 
in Fig. 13. It is clear that the northern storm contains some rich 
dynamics that cause ammonia depletion below the cloud base. 
We discuss the northern storm in more detail Section 5.2. The sub- 
tropical bands at ±9° have a lot of structure, with alternating 
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Fig. 13. (a) Map of ammonia RH in the cloud layer from March 20, 2011. Black regions would indicate supersaturation of ammonia in the cloud layer (we do not see any). 
Here, the black regions are due to the cold rings blocking the emission from the atmosphere. Green to blue regions are regions that require ammonia depletion below the 
ammonia cloud layer (i.e. RH < 0 in the cloud layer). The northern storm is blue, indicating low ammonia concentrations in the storm that extend to layers beneath the clouds. 
There are many local regions in the subtropical bands, as well as a storm in the southern hemisphere near 325°W, -43°, that require ammonia depletion below the clouds, (b) 
Two Cassini ISS images of the northern storm from Fig. 4 of Sayanagi et al. (2013), with the corresponding parts of the 2.2-cm maps show beneath each of them. These are the 
closest dates we have to the 2.2-cm map date. Cloud heights are distinguished by the three color filters - red (CB2 - 750 nm), green (MT2 - 727 nm), and blue (MT3 - 
889 nm). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) 


regions of high and low ammonia RH. In every map, there are local 
regions in the two subtropical bands that require ammonia deple- 
tion below the cloud base. As seen by the red-orange color in the 
figure, the band from -15° to -30° is more humid than other lat- 
itudes. The narrow bands at -33° and -37° are drier than their 
surrounding latitudes. South of these bands, there is a storm near 
325°W, -45° that requires ammonia depletion below the clouds. 

For comparison, we look at ammonia abundances reported by 
Fletcher et al. (2011a). Their values were derived from Cassini- 
VIMS 4.6-5. 1 pm spectra taken in April 2006. In Fig. 14i of their pa- 
per, the ammonia mole fraction is given as a function of latitude for 
their sensitivity range of 1-4 bars. There is a drastic difference 
between the equator and just off the equator in both hemispheres, 
with high ammonia abundance centered on the equator and 
extending to about ±5° planetocentric (±6° planetographic), 
and relatively low abundance at ±10° planetocentric (±12° 
planetographic). Our Fig. 9 is qualitatively consistent with these re- 
sults— there is high ammonia at the equator (low 2.2-cm bright- 
ness temperature) and low ammonia in the subtropical bands. 
However, our figure suggests that there would be larger dips in 
the ammonia abundance in the subtropical bands than is shown 
in Fletcher et al. (2011a). Their Fig. 14i shows NH 3 mole 
fractions of 450 ppm at the equator, 110 ppm at ±8-10°, and 
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Fig. 14. Saturn zonal wind velocity (solid line) and mean residual temperature for 
the March 2011 map averaged over the storm longitudes (dashed line) as a function 
of planetographic latitude. From -4° to +4° the single December 2009 RPC scan is 
used (+ signs), as in Fig. 5. The zonal velocity profile used in this figure is from 
Cassini ISS data (Garda-Melendo et al., 2011). 
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150-200 ppm at higher latitudes. Relative to our 3 x solar reference 
state and not correcting for the removal of ammonia due to the 
NH 4 SH cloud, these correspond to DF 5bar =1.3 at the equator, 
0.31 at ±8-10°, and 0.42-0.56 at high latitudes. Our data require 
DF 5 bar = 1 at the equator, 0.45-0.55 at ±8-10° and 0.8-0.9 at higher 
latitudes (also not correcting for the NH 4 SH cloud). Thus Fletcher 
et al. observe less ammonia in the subtropics and high latitudes 
and more ammonia at the equator than we do. They also see an in- 
crease in the ammonia abundance in the southern hemisphere be- 
tween -20° and -30° which is consistent with the dip in 
brightness temperature that we observe there (Fig. 9). Given the 
differences in measurement techniques and altitudes covered, 
our results tend to agree with those of Fletcher et al. (2011a). 

What dynamical mechanisms could cause a latitudinal thermal 
emission profile like the one we see in Fig. 9? We consider three 
possibilities. The first is that the brightness temperature variations 
are due to real latitudinal temperature variations, with the ammo- 
nia mixing ratio independent of latitude up to cloud base and sat- 
urated above. However, Fletcher et al. (2007) derived latitudinal 
temperature gradients on Saturn based on Cassini/CIRS observa- 
tions, and the temperature contours near the equator are nearly 
flat, at least below the 0.5-bar level (their Fig. 2). In contrast, we 
observe brightness temperature fluctuations on the order of 10I< 
between the bands at ±9° and the equator. Fletcher et al. (2007) 
do observe a temperature dip in the equatorial region above the 
0.4 bar level, with colder temperatures at the equator than in the 
subtropics by about 5 l< in the 0.2-0.4 bar region. This would be 
consistent with the brightness temperature pattern that we ob- 
serve, but at altitudes above the 0.4-bar pressure level the density 
of ammonia gas is so small that we are not sensitive to these upper 
levels. The C1RS observations only provide temperatures above the 
1 bar level so we cannot compare with the levels probed at 2.2 cm. 

In the cloud region, the ammonia concentration falls off accord- 
ing to the saturation vapor pressure (the Clausius-Clapeyron rela- 
tion). When the atmospheric temperature decreases in this region, 
holding the RH constant at saturation, the ammonia concentration 
decreases and the weighting function moves to deeper (and war- 
mer) levels, offsetting the decrease in temperature. Thus bright- 
ness temperature variations are buffered in the cloud layer, and 
we would not expect to see a large change in brightness tempera- 
ture due to a change in atmospheric temperature. This is the same 
effect that we discussed in Section 2. To confirm this, we used our 
radiative transfer model to test the sensitivity of the 2.2-cm bright- 
ness temperature to the model reference temperature at the 1-bar 
level. As we vary the reference temperature, the atmosphere re- 
mains saturated above the cloud base and the entire temperature 
profile shifts by approximately the same value as at the 1-bar level, 
including at the weighting function peak. In order to get brightness 
temperature variations on the order of 10 K, the 1-bar temperature 
needs to be varied by 50 K. The CIRS observations seem to rule out 
such large temperature swings from latitude to latitude, so we con- 
sider other options. 

The second possibility is that upwelling in the equatorial region 
and downwelling on either side in each hemisphere produces an 
ammonia vapor distribution that is compatible with our observa- 
tions. We postulate that air upwells at the equator, advecting 
ammonia-rich air from below, as suggested by previous authors 
to explain equatorial winds, composition, and clouds observed by 
Cassini (i.e. Yamakazi et al., 2005; Fletcher et al., 2011a). The 
ammonia precipitates out in the updrafts, and the dry air moves 
poleward, descending at latitudes out to ±9° in each hemisphere. 
This is the Hadley cell model. The rings of air moving poleward 
have to lose angular momentum to avoid spinning up while their 
distance to the rotation axis is decreasing. In other words, there 
has to be an eddy momentum flux (EMF) divergence within ±9° 
of the equator. On Earth, the region of EMF divergence extends to 


±30°, which is the poleward edge of the Hadley cell. The subtropics, 
which occupy the bands from 10° to 30° in each hemisphere, are 
marked by net downwelling and generally low relative humidity. 
They are “dry,” like the bands near ±9° on Saturn. On Earth the tro- 
pospheric jet streams are located in the bands from 30° to 40° in 
each hemisphere, and the eddies, which arise from instability of 
the jet streams, are responsible for the EMF divergence at lower 
latitudes. 

There are differences between the Hadley circulation on Earth 
and that on Saturn. First, in the published data there is only a hint 
of a zonal wind maximum, i.e., a jet stream, at ±10° on Saturn (Gar- 
cia-Melendo et al., 2011), as shown in Fig. 14. At the equator there 
is a zonal wind maximum, which has no terrestrial analog. 
Poleward of ±10° there is EMF convergence, which pumps the 
equatorially superrotating jets on both Saturn and Jupiter (e.g. 
Ingersoll et al., 1981; Salyk et al., 2006; Del Genio et al., 2007). 
The Hadley cell model requires that the convergence become 
divergence within ±9° of the equator. This seems to work for Jupi- 
ter, which has zonal wind maxima at ±6-7° (planetocentric) and 
EMF divergence between ±5° (planetocentric), i.e., u'v' increasing 
with latitude from -5° to +5° (see Fig. 5 of Salyk et al., 2006). 
Whether it works for Saturn is uncertain, because trackable cloud 
features are scarce close to the equator and it has been impossible 
to measure the EMF there. 

Lack of a solid surface to add angular momentum on the return 
flow is another difference between Earth and Saturn. On Earth, the 
low latitude easterlies gain westerly angular momentum from the 
surface as the air moves toward the equator. It is not clear that this 
would happen on a fluid planet. The return flow could be at any 
depth, and we do not know if rings of fluid exchange angular 
momentum there or not. Right now, the pieces of evidence for 
the Hadley cell model on Saturn are the high ammonia abundance 
at the equator, the low ammonia abundance and hint of zonal 
velocity maxima near ±9°, and the observed EMF divergence in 
the band at Jupiter’s equator. 

Fletcher et al. (2011a) present evidence that two stacked merid- 
ional circulation cells, rotating in opposite directions, exist in Sat- 
urn’s troposphere. From analysis of Cassini VIMS data, they find 
that deep PH 3 and AsH 3 show local maxima on either side of the 
equator, whereas the PH 3 scale height, the upper cloud opacity, 
and the NH 3 show local maxima at the equator. Stacked cells seem 
to explain the contradictory evidence of upwelling and downwel- 
ling both on and off the equator. As Fletcher et al. point out, the 
stacked cell hypothesis was originally invoked at Jupiter (Ingersoll 
et al„ 2000; Gierasch et al., 2000; Showman and de Pater, 2005), 
where upper cloud opacity and NH 3 indicated upwelling in the 
zones, but lightning and other evidence of moist convection indi- 
cated upwelling in the belts. Other authors also explore the possi- 
bility of meridionally overturning cells in Saturn’s atmosphere (i.e. 
Del Genio et al., 2009). Since the present paper mainly concerns 
NH 3 , we do not attempt to synthesize all the Saturn data at this 
time. 

5.2. 2030-2033 Northern storm 

The third possibility is that there is some process in giant planet 
atmospheres that causes the “drying out” following convective 
events. The northern storm in the March 201 1 map is very dry with 
respect to ammonia vapor (Fig. 13). What is it about the storm that 
produces such low ammonia RH? The head of the storm is at 40°, 
180°W (Sayanagi et al., 2013), and is indicated by the red triangle 
at the top of the figure in Fig. 13b. The green to blue trail to the 
west of the head is the tail of the storm, which wrapped all the 
way around the planet until it collided with the head in February 
2011 (Fischer et al., 2011). By March 20, 2011, when this 2.2-cm 
data was collected, the storm had been in existence for almost 
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4 months. The storm was a copious producer of lightning, which 
was detected at radio frequencies (Fischer et al., 2011) and in vis- 
ible light (Dyudina et al., in preparation). Fig. 13 shows that, 
although the storm is very dry everywhere (RH ^ 0, Fig. 13a), there 
is also quite a bit of brightness temperature structure within the 
storm, which may be compatible with the Sanchez-Lavega et al. 
(2012) description of the three branches of the storm and the 
wave-like patterns within the storm tail. 

Fig. 13b shows images taken by the Cassini imaging science 
subsystem (ISS) on March 7 and March 17, 2011 (Sayanagi et al., 
2013), with the corresponding pieces of the RADAR maps below 
them (from March 20, 2011). The top panel of Fig. 13b spans 
0-200°W longitude, and the bottom panel spans 60-1 30°W longi- 
tude. Both panels show 24.5-45° planetographic latitude (20.5-40° 
planetocentric). Red, green and blue color channels correspond to 
Cassini ISS camera’s CB2 (750 nm), MT2 (727 nm), and MT3 
(889 nm) filters, respectively, and they convey the cloud top 
heights. The altitude of the features generally increases in the order 
of red-green-blue (Sayanagi et al., 2013). There are no images clo- 
ser in time to the March 20, 2011 radiometer image in Fig. 13a, for 
several reasons. Most important, the two instruments point along 
different axes of the spacecraft and cannot take data simulta- 
neously. Also, the 2.2-cm map is taken near closest approach, when 
the field of view of the wide-angle camera is 15-20° latitude, 
which is a small fraction of the planet. Finally, the observations 
were scheduled months in advance, long before the advent of the 
storm. Thus the detailed features in Fig. 13a are not necessarily 
the same features that appear in the ISS images in Fig. 13b. Some 
features, however, seem to be captured by the RADAR map, for 
example the anticyclonic vortex appears as a circular region of 
low RH (dark blue) near 15°W. For an example of the changes that 
take place over an 11-h period, see Fischer et al. (2011). 

Fig. 14 shows the zonal wind velocity (solid line) as a function 
of latitude with the March 2011 residual temperature averaged 
over the longitudes of the storm (dashed line) overlaid on it. The + - 
signs between -4° and 4° depict the December 2009 RPC scan 
brightness residuals used in Fig. 9. Storm alley in the southern 
hemisphere lies in the westward jet between -40° and -45°. The 
northern storm is also located in a westward jet, but in the oppo- 
site hemisphere. The head of the storm is centered in the middle of 
the westward jet near 41 ° (35° planetocentric, Fischer et al., 201 1 ). 
Some of the depletion of ammonia vapor by the northern storm 
gets pushed slightly to the north of the westward jet, which is con- 
sistent with work by Fletcher et al. (2011b) and Sanchez-Lavega 
et al. (2012), who observe two eastward branches to the north 
and south of the storm, with ammonia vapor depletion in the 
southern branch (Fletcher et al., 2011b). This pattern may be anal- 
ogous to the storms in the southern hemisphere which are to the 
south of the westward jet (Porco et al., 2005). We suspect that 
the mechanism for “drying out” the small storms in the southern 
hemisphere and the 2010-2011 northern storm are the same. 
Ammonia depletion is consistent with the ISS (imaging science 
subsystem) team’s interpretation of their observations of the 
southern storms at near-IR wavelengths. To measure clouds at var- 
ious altitudes, they use three filters, like in Fig. 13b, where the 
absorption by methane gas is strong, weak, and negligible, respec- 
tively (Porco et al., 2005; Dyudina et al., 2007). From their observa- 
tions they infer that on the first day the storms are active 
convective regions of optically thick, high clouds. On the third or 
fourth day they become stable circular clouds with optically thin, 
high haze above dark regions with no deep clouds - they become 
holes in the clouds. If the clouds are made of ammonia ice, then 
a hole in the clouds is consistent with depletion of ammonia vapor. 

The VIMS (visible and infrared imaging spectrometer) team 
uses 352-bandpass spectra ranging from 0.35 to 5.1 pm, although 
the instrument has other ranges and resolutions. The dark regions 


seen in the near-IR are dark at all wavelengths, and Baines et al. 
(2009) interpret them as carbon-impregnated water frost rather 
than absence of deep clouds, the latter being the ISS interpretation. 
The carbon could come from dissociation of methane by lightning, 
since these are lightning clouds (Dyudina et al., 2007; Fischer et al., 
2011 ). 

Both the ISS description and the VIMS description seem to ex- 
plain the dark color of the convective clouds after the third day. 
The low RH of ammonia inferred from the 2.2-cm observations fa- 
vors a hole in the clouds if the cloud particles are ammonia ice. The 
2010-2011 northern storm, whose head resembles the high, thick 
convective clouds in ISS images, and whose tail is depleted in 
ammonia in the 2.2-cm maps, also supports the hole-in-the-cloud 
interpretation. Of course, the tail could be depleted in ammonia 
clouds and ammonia vapor but still have carbon-impregnated 
water ice clouds. The VIMS spectra of the head and tail of the 
northern storm will help resolve these different interpretations. 

Convective storms seem to evolve into ammonia-poor regions, 
both in the southern hemisphere hot spots and in the northern 
storm. Also, the remnant of the northern storm is bright at 5 pm, 
which signals “an unusual dearth of deep clouds” (Baines, private 
communication, 2012). This agrees with the ISS interpretation of 
the southern hemisphere dark spots - that deep clouds are absent. 
So the question arises, why should convection “dry out" the atmo- 
sphere - removing gaseous ammonia and deep clouds? The answer 
may be in the nature of convection in hydrogen-helium atmo- 
spheres, in particular in the effect of mass loading by molecules 
of NH 3 , H 2 S, and H 2 0, which have higher molecular mass than 
the ambient atmosphere. The parcels that have lost their load of 
the more massive molecules will have the most buoyancy, and 
they might rise the highest. These parcels would look “dry” when 
viewed from the top of the atmosphere. Mass loading might also 
explain the intermittency of convection on Saturn, but the details 
have yet to be worked out. 

There are several models of moist convection in giant planet 
atmospheres. Some are axially symmetric (Yair et al., 1995a,b), 
and some are three dimensional (Hueso and Sanchez-Lavega, 
2004), with the possibility of wind shear and precipitation on 
one side of the central updraft. These models start with an unstable 
initial state and follow the convective plume as it develops over a 
period of several hours or 1 day. Other models incorporate 
cumulus parameterizations into giant planet general circulation 
models (Del Genio and McGrattan, 1990; Palotai and Dowling, 
2008). Sanchez-Lavega et al. (2012) studied a mass source in a 
shear flow patterned after the westward jet at 40° planetographic 
that produces a GWS (Great White Spot) long tail. To the best of 
our knowledge, none of these models capture the 20-30 year cycle 
of planet-encircling storms or the ammonia depletion following 
convective events. 


6. Conclusions 

J13 used the 2.2-cm brightness temperature observations of 
Saturn presented in conjunction with radiative transfer calcula- 
tions to produce residual brightness temperature maps for five 
dates between 2005 and 2011. In this work we analyzed the maps 
by making adjustments to the vertical ammonia distribution using 
the JAMRT program. We find that ammonia vapor must be de- 
pleted below the cloud base in some regions in order to obtain 
temperatures in agreement with observations. The observed 
brightness temperatures are consistent with a deep abundance of 
3-4x solar ammonia (3. 6-4.8 x 10~ 4 volume mixing ratio) with 
varying depletion factors relative to the standard model, which is 
saturated above cloud base. The depletion must extend to 2 bars 
or deeper for brightness temperatures >160 K. To obtain these 
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results, we assume that Saturn’s latitudinal temperature profile is 
constant in our sensitivity range of 0.5-2 bars. The highest bright- 
ness temperatures we see are in the 2010-2011 northern storm 
(1 65.7 K) and in the subtropical latitudes of the October 2009 
map (167 K). The most striking feature, evident in Fig. 9, is the dif- 
ference in the 2.2-cm brightness temperature right at the equator 
versus that just off the equator. This implies that there are some 
interesting atmospheric dynamics at play in the equatorial region. 

We presented three options to explain the brightness tempera- 
ture pattern observed at 2.2-cm. One is that it represents real tem- 
perature variations, i.e., temperature variations from place to place 
at constant levels in the atmosphere. In principle, brightness 
temperature variations can be due to both physical temperature 
and absorber concentrations; however, the buffering effect makes 
the former option unlikely in our case. The second option is that 
large-scale upwelling and downwelling, like the Earth’s Hadley cir- 
culation, creates dry zones like the subtropics of Earth. However 
this option does not explain the drying that follows the small-scale 
convective events in the southern hemisphere. The third option, 
which we do not explore in any depth, is that the drying is an 
intrinsic property of convection in giant planet atmospheres, and 
that it applies not just to the small southern lightning storms but 
also to the northern storm of 2010-2011. The subtropical dry 
bands are not copious producers of lightning and moist convection, 
so the third option might not apply there. It may that that the 
meridional circulation explanation applies to the subtropics, and 
the deep convection explanation applies to the lightning storms. 

Modeling of atmospheric circulations on giant planets has pro- 
ven difficult, and the lack of observational data in the deep atmo- 
sphere below the clouds is a limitation. It is important to try to 
reconcile the atmospheric motions with latitudinal distributions 
of tropospheric gases such as NH 3 and PH 3 . This paper provides 
some insight into the distribution of ammonia vapor, with the 
hope that more work can be done with these observations to rec- 
oncile them with the energy and momentum balances of Saturn’s 
atmosphere. 
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